Journal of The Korean Astronomical Society 
37: 1 ~ 5, 2004 



The quest for cosmic ray protons in galaxy clusters 

C. Pfrommer and T. A. Ensslin 
Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1 , PO Box 1317, 85741 Garching, Germany 
E-mail: pfrommer@mpa-garching.mpg.de, ensslin@mpa-garching.mpg.de 
(Received October 31, 2004; Accepted December 1, 2004) 



ABSTRACT 



o 
o 

O 

Q 

in 



m 

(N 

o 

Of 

6 

c3 



There have been many speculations about the presence of cosmic ray protons (CRps) in galaxy clusters over the past two 
decades. However, no direct evidence such as the characteristic y-ray signature of decaying pions has been found so far. 
These pions would be a direct tracer of hadronic CRp interactions with the ambient thermal gas also yielding observable 
synchrotron and inverse Compton emission by additionally produced secondary electrons. The obvious question concerns the 
type of galaxy clusters most likely to yield a signal: Particularly suited sites should be cluster cooling cores due to their high 
gas and magnetic energy densities. We studied a nearby sample of clusters evincing cooling cores in order to place stringent 
limits on the cluster CRp population by using non-detections of EGRET. In this context, we examined the possibility of a 
hadronic origin of Coma-sized radio halos as well as radio mini-halos. Especially for mini-halos, strong clues are provided 
by the very plausible small amount of required CRp energy density and a matching radio profile. Introducing the hadronic 
minimum energy criterion, we show that the energetically favored CRp energy density is constrained to 2% + 1% of the 
thermal energy density in Perseus. We also studied the CRp population within the cooling core region of Virgo using the 
TeV y-ray detection of M 87 by HEGRA. Both the expected radial y-ray profile and the required amount of CRp support this 
hadronic scenario. 
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I. Galaxy clusters as a window to non-thermal intra- 
cluster medium components 

Cooling cores in galaxy clusters are especially well 
suited places to find traces of otherwise nearly invisible non- 
thermal components of the intra-cluster medium (ICM) due 
to the extreme gas densities observed in these central re- 
gions. The faded and therefore invisible remnants of ra- 
dio galaxy cocoons, so-called radio ghosts (EnBlin 1999) or 
ghost cavities, were first detected in cooling cores by the 
absence of X-ray emissivity in the ghost's volume in con- 
trast to the highly X-ray luminous cooling cores gas sur- 
rounding it (Bohringer et al. 1993; Fabian et al. 2000, and 
many recent Chandra observations). Cosmic ray electrons 
(CRes) are seen in cooling cores by their radio synchrotron 
radiation within strong magnetic fields. Cosmic ray protons 
(CRps) in the ICM are most likely to be detected for the first 
time within cooling cores via their hadronic interaction with 
the dense cooling core gas leading to y-rays and CRes. 

A better knowledge of these non-thermal components of 
the ICM - especially in the cooling core regions - is highly 
desirable, since they play important roles in the heat balance 
of the gas through heating by CRps, radio ghost buoyant 
movements, and suppression of heat conduction by mag- 
netic fields. Additionally, such non-thermal components are 
tracers of the violent dynamics of the ICM and may help to 
solve some of the puzzles about cooling cores. 

In this article, we present our recent progress in con- 
straining CRps by using limits on the y-ray emission as 
well as the observed ra dio e mission of nearby galaxy clus- 
ters (Sect. HJl] and Sect. IIVI . Finally, we present the novel 
concept of the hadronic minimum energy condition which 



should serve as unambiguous energetic criterion for the ap- 
plicability of the hadronic model of radio synchrotron emis- 
sion in galaxy clusters (Sect.|VJ. We focus on ideas and re- 
sults, leaving the technical details to the publications given 
in the reference list. 

II. Hadronic interactions of cosmic ray protons 

Approximately once in a Hubble time, a CRp collides 
inelastically with a nucleon of the ambient ICM gas of non- 
cooling core clusters. Within cooling cores, such collisions 
are much more frequent due to the higher target densities. 
Such inelastic proton (p) nucleon (AO collisions hadroni- 
cally produce secondary particles like relativistic electrons, 
positrons, neutrinos and y-rays according to the following 
reaction chain: 



p+N 

_± 



2N + n ±/0 
2y. 



+ V e /V e + Vn + V„ 



The resulting y-rays can be detected directly with current 
and future y-ray telescopes. The relativistic electrons and 
positrons (summarized as CRes) are visible due to two radi- 
ation processes: inverse Compton scattering of background 
photon fields (mainly the cosmic microwave background, 
but also starlight photons) and radio synchrotron emission in 
ICM magnetic fields. Especially the latter process provides 
a very sensitive observational signature of the presence of 
CRps in cooling cores not only because of the tremendous 
collecting area of radio telescopes, but also due to the strong 
magnetic fields in cooling cores, as Faraday rotation mea- 
surements demonstrate. 
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Fig. 1. — Theoretical y-ray spectra of the Perseus cluster 
resulting from CRp hadronic interactions with the cooling 
core gas. The assumed CRp spectra are normalized to be 
marginally consistent with the EGRET E y > 100 MeV non- 
detection limits. Results for different proton spectral indices 
are shown. The 7r°-decay bump is clearly visible. Secondary 
CRes scatter cosmic microwave background photons into 
the y-ray regime by the inverse Compton process. The re- 
sulting power-law emission dominates the low energy part 
of the spectrum. The CRe population was calculated ne- 
glecting synchrotron radiation energy losses. Thus, realistic 
inverse Compton spectra due to this process will exhibit a 
lower normalization than displayed here. 

III. Gamma-ray emission 

(a) Gamma-ray constraints from EGRET 

Assuming that the CRp population can be described by a 
power law distribution in momentum, Pfrommer & EnBlin 
(2004a) developed an analytical formalism to describe the 
secondary emission spectra from hadronic CRp interactions 
which exhibit the simplicity of textbook formulae. This for- 
malism was applied to a sample of nearby X-ray luminous 
galaxy clusters which are also believed to be powerful y-ray 
emitters owing to the present high target densities. Synthetic 
y-ray spectra of the Perseus cooling core cluster, which are 
calculated using this formalism, are shown in Fig. ^ As- 
suming that the CRp population follows the spatial distribu- 
tion of the thermal ambient intra-cluster gas, we can define 
the CRp scaling ratio 



where ecRp and e t h denote the CRp and thermal energy den- 
sity, respectively. This scaling effectively implies an iso- 
baric compression of the CRps during the formation of the 
cooling core. The parent CRp spectra, which give rise to the 
hadronically induced y-ray spectra, are normalized to upper 
limits on the y-ray emission obtained by EGRET observa- 
tions for energies E 7 > 100 MeV (Reimer et al. 2003). For 
nearby cooling core clusters, this analysis constrains CRp 



Fig. 2. — EGRET y-ray constraints on the (central) CRp 
energy density (ecRp, for various CRp spectral indices) in 
terms of the cluster central thermal energy density (eth) 
for cooling core and non-cooling core clusters. The con- 
straint from the Perseus cluster radio mini-halo observation 
is also displayed while assuming a typical cooling core clus- 
ter magnetic field strength (blue arrow). If radio mini-halos 
are entirely produced hadronically by CRps then the con- 
straint derived for the Perseus cluster is an actual measure- 
ment for XcRp and not only an upper limit. 

energy densities relative to the thermal energy density to 
^CRp ~ 10%. For the full sample of nearby X-ray lumi- 
nous clusters, upper limits of the same order of magnitude 
were obtained (cf. Fig.|5J. It is obvious from this compila- 
tion that cooling core clusters are extremely well suited to 
visualize even small CRp populations. 

The real y-ray spectra are not expected to be far below 
these spectra: many processes like supernova driven galac- 
tic winds, structure formation shock waves, radio galaxy ac- 
tivity, and in-situ turbulent particle acceleration support this 
expectation. These processes should have produced a CRp 
population characterized by an energy density relative to the 
thermal energy density of at least a few percent. 

(b) Gamma-rays from the central region of the Virgo 
cluster 

Recently, the HEGRA collaboration (Aharonian et al. 
2003) announced a TeV y-ray detection from the giant el- 
liptical galaxy M 87 which is situated at the center of the 
cooling core region of the Virgo cluster. Using imaging at- 
mospheric Cerenkov techniques, this y-ray detection was 
obtained at a 4-cr significance level. On the basis of their 
limited event statistics, it is inconclusive whether the de- 
tected emission originates from a point source or an ex- 
tended source. Despite testing for intermittent behavior of 
M 87, no time variation of the TeV y-ray flux has been 
found. 

Pfrommer & EnBlin (2003) ap plied the previously de- 
scribed analytical formalism (Sect. |(a)) of secondary y ray 
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Fig. 3. — Upper panel: Modeled y-ray surface flux pro- 
file of the cooling core region within the Virgo cluster. The 
profile is normalized by comparing the integrated y-ray flux 
above 730 GeV to HEGRA data of M 87 within the inner- 
most two data points. Lower panel: Comparison of de- 
tected to predicted y-ray flux within the central aperture and 
the innermost annuli for two different widths of the PSF. 
The red lines correspond to an optimistic PSF of cr — 0.05° 
whereas the green lines are calculated for a conservative 
cr = 0.08° for a soft y-ray spectrum. The vertical dashed 
line separates the data from the noise level at a position cor- 
responding to r ± — 37.5 kpc. 



emission spectra resulting from hadronic CRp interactions 
to the central cooling core region of the Virgo cluster. While 
combining the observed TeV y-ray emission to EGRET up- 
per limits, they constrain the CRp spectral index a<?y < 2.3 
provided the y-ray emission is of hadronic origin and the 
population is described by a single power-law ranging from 
the GeV to TeV energy regime. 

A comparison of the observed to the predicted y-ray 
emissivity profiles is shown in Fig. [3] In order to allow for 
finite resolution of the Cerenkov telescope, the real profile 
was convolved with the point spread function (PSF) of the 



Fig. 4. — The Perseus radio mini-halo at 1 .4 GHz from Ped- 
lar et al. (1990) with an extent of 160 /i^J kpc in diameter. 
At the very center is an extremely bright flat-spectrum core 
owing to relativistic outflows of the radio galaxy NGC 1275. 
This outflow was subtracted from this image to make the 
extended structure more visible. The position of the radio 
galaxy is marked by the green cross in the center of the red 
region. 



HEGRA instrument. The more optimistic PSF of cr — 0.05° 
for harder y-ray spectra being favored by our model as well 
as the conservative choice of a PSF derived from softer 
Crab-like spectra are both consistent with the observed data. 
Considering an aged CRp population which has already suf- 
fered from significant Coulomb losses, CRp scaling ratios 
of the order of Xcr p ~ 50% are obtained for the innermost 
region of Virgo within r ± = 37.5 kpc. 

Since the emission region is dominated by the giant radio 
galaxy M 87, other mechanisms like processed radiation of 
the relativistic outflow or dark matter annihilation could also 
give rise to the observed y-ray emission. Nevertheless, the 
hadronic scenario probes the CRp population within a mix- 
ture of the inter stellar medium of M 87 and the ICM of the 
center of Virgo yielding either an upper limit or a detection 
on the CRp population, provided this scenario applies. 

IV. Radio mini-halos 

The Perseus radio mini-halo at 1.4 GHz is shown in 
Fig. as observed by Pedlar et al. (1990). The emission 
due to the relativistic jet has been subtracted to make the 
extended structure more visible. The spatial extent of the 
radio mini-halo of 160 hjl kpc in diameter is too large to be 
accounted for by synchrotron emission of direct accelerated 
CRes in structure formation or accretion shocks. Thus one 
needs to consider other injection processes of CRes into the 
ICM which are responsible for the observed diffuse radio 
emission. Besides the reacceleration scenario of mildly rel- 
ativistic CRes (y 100 - 300) which are accelerated in-situ 
by turbulent Alven waves (Gitti et al. 2002), the hadronic 
injection scenario of CRes is a promising alternative. 
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Fig. 5. — Radio emissivity profile of the Perseus mini-halo 
at 1.4 GHz observed by Pedlar et al. (1990) (crosses) and 
predicted due to hadronic interactions of CRp by Pfrom- 
mer & EnBlin (2004a). Both the required CRp scaling ra- 
tio XcRp - 2% and the morphological match of the ob- 
served and predicted radio emissivities strongly indicate the 
hadronic origin of radio mini-halos. 



Azimuthally averaging the radio emission shown in Fig.|4] 
yields the radio emissivity profile as displayed in Fig. [5] 
For comparison, we overlaid synthetic radio surface bright- 
ness profiles resulting from the hadronic scenario. There is 
a perfect morphological match between observed and pre- 
dicted profiles. Assuming a typical magnetic field strength 
for cooling cores of lOyi/G, we derive an upper limit for the 
CRp population of Xqr p - 2% by normalizing our simu- 
lated to the observed profile (cf. Fig. |2j. 

The radio emissivity of hadronically generated CRes de- 
pends on the assumed magnetic field strength. Thus, upper 
limits on Xcr p rely on the same assumption. However, in 
the case of strong magnetic fields (above 3 /iG), the depen- 
dence of XcRp on the assumed field strength is very weak, 
since synchrotron losses dominate in this regime. If ra- 
dio mini-halos are entirely produced hadronically by CRps, 
then the constraint derived for the Perseus cluster is an ac- 
tual measurement for Xcr p and not only an upper limit. 
Alongside the perfect morphological match of the radio syn- 
chrotron profiles, the comparably small number of CRps re- 
quired to account for the observed radio mini-halo suggests 
a hadronic origin of the radio mini-halo in Perseus. 

V. Hadronic minimum energy criteria 

As an extension of the previous work, we estimated mag- 
netic field strengths of radio emitting galaxy clusters by 
minimizing the non-thermal energy density — contained 
in relativistic electrons, protons, and magnetic fields — 
with respect to the magnetic field strength (Pfrommer & 
EnBlin 2004b). As one boundary condition, the implied syn- 
chrotron emissivity is required to match the observed value. 
Additionally, a second boundary condition is required math- 



ematically which couples CRps and CRes. For the classical 
case, a constant scaling factor between CRp and CRe en- 
ergy densities is assumed. However, if the physical connec- 
tion between CRps and CRes is known or assumed, a phys- 
ically better motivated criterion can be formulated. As such 
a case, we introduce the minimum energy criterion within 
the scenario of hadronically generated CRes. The work by 
Pfrommer & EnBlin (2004b) provides simple self-consistent 
recipes for applying the classical and hadronic minimum en- 
ergy criterion in typical observational situations. 

Alongside, we provide theoretically expected tolerance 
regions which measure the deviation from the minimum en- 
ergy states by one e-fold: We use logarithmic measures of 
the curvature radius at the extremal values in order to char- 
acterize the 'sharpness' of the minima. These regions have 
the meaning of a quasi-optimal realization of the particular 
energy densities. 

The philosophy of this approach is to provide a crite- 
rion for the energetically least expensive radio synchrotron 
emission model possible for a given physically motivated 
scenario. There is no first principle enforcing this mini- 
mum state to be realized in nature. However, our mini- 
mum energy estimates are interesting in two respects: First, 
these estimates allow scrutinizing the hadronic model for 
extended radio synchrotron emission in clusters of galax- 
ies. If it turns out that the required minimum non-thermal 
energy densities are too large compared to the thermal en- 
ergy density, the hadronic scenario will become implausible 
to account for the extended diffuse radio emission. In this 
respect, our criteria is a way to test the robustness of the 
previous results with respect to the observationally available 
parameter space spanned by the CRp spectral index and the 
(unknown) distribution of the magnetic field strength. Sec- 
ondly, should the hadronic scenario be confirmed, the min- 
imum energy estimates allow testing for the realization of 
the minimum energy state for a given independent measure- 
ment of the magnetic field strength. 

Application to the radio halo of the Coma cluster (Deiss 
et al. 1997) and the radio mini-halo of the Perseus clus- 
ter (Pedlar et al. 1990) yields equipartition between cosmic 
rays and magnetic fields within the expected tolerance re- 
gions. In the hadronic scenario, the inferred central mag- 
netic field strength ranges from 2.4 /iG (Coma) to 8.8 fiG 
(Perseus), while the optimal CRp energy density is con- 
strained to 2% + 1% of the thermal energy density (Perseus) 
(cf. Fig.|6j. Pfrommer & EnBlin (2004b) discuss the possi- 
bility of a hadronic origin of the Coma radio halo while cur- 
rent observations favor such a scenario for the Perseus radio 
mini-halo. Combining future expected detections of radio 
synchrotron, hard X-ray inverse Compton, and hadronically 
induced y-ray emission should allow an estimate of volume 
averaged cluster magnetic fields and provide information 
about their dynamical state. 

VI. Conclusions 

• We argue that cooling cores of galaxy clusters are 
well suited to reveal or constrain any cosmic ray pro- 
ton population via radiation from hadronic interac- 
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Fig. 6. — Profiles of the CRp-to-thermal energy density XcR Pmia (r) (solid) and magnetic-to-thermal energy density X Bmin (r) 
(dotted) as a function of deprojected radius are shown. The different energy densities are obtained by means of the hadronic 
minimum energy criterion. The left panel shows profiles of the Coma cluster while the right panel represents profiles of 
the Perseus cluster. The light shaded areas represent the logarithmic tolerance regions of XB m ,„(r) and XcR Pmin (r), respectively, 
while the dark shaded regions indicate the overlap and thus the possible equipartition regions in the quasi-optimal case. 



tions with the ambient gas nuclei. Such collisions 
lead to y-rays and cosmic ray electrons. The former 
would have been seen above 100 MeV by the EGRET 
telescope if the cosmic ray protons had energy densi- 
ties relative to the thermal gas exceeding 10% (Pfrom- 
mer & Enfilin 2004a). 

• The giant elliptical galaxy M 87 in the center of the 
Virgo cluster cooling core region has recently been 
detected at TeV energies by the HEGRA instrument 
(Aharonian et al. 2003). These y-rays could be pro- 
duced by hadronic interactions of a cosmic ray proton 
population if its spectral index a<?y < 2.3 and its en- 
ergy density is of the order of 50% of the gas within 
the transition/mixture of inter-stellar and intra-cluster 
medium within M 87 (Pfrommer & Enfilin 2003). 

• Cosmic ray electrons produced in hadronic interac- 
tions are a very sensitive indicator of cosmic ray 
protons due to their strong emissivity. Radio syn- 
chrotron emission of such electrons in strong cooling 
core magnetic fields of ~ 10/iG limit the cosmic ray 
proton energy density to ~ 2% or less compared to 
the thermal one in the Perseus cluster cooling core 
(Pfrommer & Enfilin 2004a). 

• Diffuse radio emission from the Perseus cluster cool- 
ing core was detected - the so called Perseus radio 
mini-halo. This radio synchrotron emission may be 
induced by cosmic ray proton interactions in the intra- 
cluster medium. This scenario is strongly supported 
by the very moderate energy requirements (2% of the 
thermal energy) and the excellent agreement between 
the observed and the theoretically predicted radio sur- 
face profile (Pfrommer & Enfilin 2004a). 



• Introducing the hadronic minimum energy criterion, 
we show that the energetically favored cosmic ray 
proton energy density is constrained to 2% + 1% of 
the thermal energy density in Perseus. Application to 
the radio halo of the Coma cluster and the radio mini- 
halo of the Perseus cluster yields equipartition be- 
tween cosmic rays and magnetic fields within the ex- 
pected tolerance regions (Pfrommer & Enfilin 2004b). 
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